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I solutions expressed in terms of Riemann invariants. We construct an abelian distribution 

\Q . of vector fields which are symmetries of the original system of PDEs subjected to certain 

first order differential constraints. The usefulness of our approach is demonstrated on simple 
and double wave solutions of MHD equations. The paper also contains a comparison of the 
conditional symmetry method with the generalized method of characteristics. 

1 Introduction 

The objective of this paper is a development of the version of the conditional symmetry 
method (CSM) (proposed in p],[2]) for the purpose of constructing simple Riemann wave 
solutions and their superpositions (i.e. multiples waves) admitted by the equations of 
magnetohydrodynamics (MHD) in (3 + 1) dimensions. The flow under consideration is 
assumed to be ideal, nonstationary and isentropic for a compressible conductive fluid 
placed in magnetic field H. The electrical conductivity of the fluid is assumed to be 
infinitely large (i.e. the fluid is an ideal conductor o \— > oo). We restrict our analysis 
to the case in which the dissipative effects, like viscosity and thermal conductivity, are 
negligible and no external forces are considered. Under the above assumptions the MHD 
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model is governed by the system of equations 
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where we use the following notation: p and p represent the density and the pressure of 
the fluid, respectively; v = (u,v,w) and H = [H^ H 2 , H 3 ) represent the velocity of the 
fluid and the magnetic field, respectively and k is an adiabatic exponent. Without loss 
of generality in this model we can set the coefficient of the magnetic permeability to 
unity jj, e = 1. The independent variables are denoted by = (t,x,y,z) E E C 1R 4 , 
\i — 0, 1, 2, 3. The system of MHD equations fll.ip - fll.4l) is of evolutionary form. It is well 
known [3] that, if the magnetic field H obeys equation (11.51) at time t — 0, then, by virtue 
of the MHD equations fll.ll) - fll.4p . it will retain this property for all t > 0. The system 
of MHD equations fll.ll) - fll.5p is composed of nine equations involving eight dependent 
variables u = (p,p,v,H) E U C R 8 . In (3 + 1) dimensions the system of equations 
fll.ip - fll.4l) can be written in the matrix equivalent evolutionary form 



u t + J2 Ai ( u ) = > 

i=l 

where 8 by 8 matrix functions A 1 , A 2 and A 3 take the form 
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The MHD system f ll.ll) -( jl~5"]) is an hyperbolic one and it is invariant under the Galilean- 
similitude Lie algebra [I]. 

In this paper we seek for solutions describing the propagation and nonlinear superpo- 
sitions of waves which can be realized in the MHD system (1 1 . 1 j) - ( TT75T) . A wave vector of 
system ( 11.61) is a nonzero vector function 

A(u) = (A (u), ... ,Ap(u)), (1.10) 

such that 

ker (Ao (u)I + Ai (u)A l ) ^ (1.11) 

holds. This means that there exists an eight-component vector 7 = (7 1 , . . . , 7 s ) for which 
the condition 

(A (u)I + K (u)O 7 = (1.12) 

is satisfied. Note that throughout this paper we use the summation convention over the 
repeated lower and upper indices, except in cases when the index is taken in brackets. The 
necessary and sufficient condition for the existence of a nonzero solution 7 of equations 

(H32H is 

rank (A (u)I + Ai (u)^) < 8 . (1.13) 

The relations (j!.12p and fll . 13j) are called the wave relation and the dispersion relation, 
respectively [5]. The wave vector A can be written in the form A = (A , A), where 

A = (Ai, A2, A3) denotes the direction of wave propagation and A Q is the phase velocity of 
the considered wave. The dispersion relation (11 . 13j) for MHD equations (11 . ip — fll . 51) takes 
the form 



5 2 |A| 2 U 2 IA| 2 



(H-Xf 
P 



5 4 |A| 4 -5 2 |AI 2 



\H\ 
P 



-a 2 ]+a 2 ^-^ 
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0, (1.14) 



where d>|A| = A c + v ■ A denotes the wave velocity with respect to a moving fluid and 

1 

a = (up/ p) 2 is the velocity of sound. Solving the dispersion equation (jl.l4|) with respect 
to 8\X\ we obtain the following eigenf unctions 
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where e = ±1 implies that the wave propagates in the right or in the left direction 
with respect to the medium. The eigenvalues (I1.15j) - (jl.l8[) correspond to entropic waves 
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5 E , Alfven waves 5a, and magnetoacoustic slow 5s and fast 5f waves, respectively. The 
eigenvectors 7 and A, corresponding to the four types of eigenvalues fll.15p - fll.18p admitted 
by the MHD system of equations fll.ip - fll.5l) . have the following form [5], [6] 



(1) for the entropic waves, we have three types of eigenvectors 

Ei ■ i El =(l P ,-(H-\),l,h), X El = (-v-A,A), (1.19) 

E 2 : lE2 = ( 7p ,-(jy-K),7,K), X l = (-v-(cTxtf), ^xff), z = l,2 (1.20) 
where a % are two linearly independent vectors, 

E 3 : 7 E 3 = (7p, 0,0,0) , A*= (-v-e^ei), 2 = 1,2,3 (1.21) 

where e*j are three linearly independent unit vectors, 



(2) for the Alfven waves 



A- : -,= (0.0. — .//). A A = |e^-^-v-A,A], e = ±1 , (1.22) 
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(3) for the magnetoacoustic waves S and F 
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where we have used the following notation for eigenvectors 7 = (^7 P ,7 P ,7, h j. They are 

related to unknown functions u = (p, p, v, H^j G U . Similarly the wave vector A is related 
to the independent variables x = (t, x, y, z) = (t, x) G E. 



The functions 



r s (x, u) = A: (u)t + A*(u)x l , s = l, ...,8 (1.24) 



are called the Riemann invariants of the wave vectors X s (see e.g. |3J,|7J). So according 
to equation (11.241) . the MHD model fll.ll) - fll.5p admits the following Riemann invariants 
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associated with wave vectors A s 

r El (x,u) = A(u) -x- (A(u) ■ v)t , (1.25) 

r E2 (x,u) = x + y — (u + v)t, (1.26) 

r Es (x, u) = x + y + z - (u + v + w)t , (1-27) 

r A (x,u) = A (u) ■ x + | e ^— ' *j U ^ - v • A (u) | t , e = ±1 , (1.28) 

V VP J 

r s (x,u) = X(u) -x+ (S s \ A(u)| - v- A(u))t, (1.29) 

r F (x,u) = A(u) -x+ (6 F \ A(u)| -v- A(u))t. (1.30) 

For any function, / : R — >■ R 8 , the equation 

u = /(r(x,u)) (1.31) 
defines a unique function u (x) on a neighborhood of x = and the Jacobian matrix is 

du a 

^ r = 0(x)- 1 A^(u(x)) 7 f s) , ^ = 0,1,2,3; a, s = 1, ... ,8 (1.32) 
where we have introduced the notation 

0(x) = 1 - J^(x,u(x)) 7f s) (r (x,u(x))) , (1.33) 

and 

df a 

Is = J~ s ■ (1-34) 

dr s 



Note that the Jacobian matrix fll.32p has rank at most equal to one. The rank-one 
solutions of the hyperbolic system (11.61) are called simple waves and always exist [3] . This 
type of solution was introduced by S.D. Poisson [8] at the beginning of the 19th century 
in connection with the equations describing a compressible isothermal gas flow. 

In this paper we construct several classes of exact solutions in the form (11.311) for the 
MHD equations (ll.lj) - (ll.5p . In particular we focus on constructing simple wave solutions, 
scattering and nonscattering double wave solutions (rank-2 solutions) obtained by the 
conditional symmetry method described below. 

This paper is organized as follows. Section 2 presents an adapted version of the CSM. 
We use it to obtain simple waves of MHD system of equations. Section 3 contains a 
detailed description of how to construct certain classes of double wave solutions admitted 
by the MHD equations. Section 4 summarizes the obtained results and contains a 
comparison of these results with the generalized method of characteristics (GMC). 

2 Conditional symmetries and simple wave solutions of 
MHD equations. 

The methodological approach adopted here is a new variant of the CSM presented in pQ 
and [2j. The notion of conditional symmetries evolved in the process of extending the 
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classical Lie theory of symmetries to partial differential equations (PDEs). This approach 
consists basically in modifying the original system by adding to it certain differential 
constraints of the first order for which a symmetry criterion is identically satisfied [9]. 
The overdetermined system of equations obtained in this way admits, in some cases, 
a larger class of Lie point symmetries and, consequently, it can provide new classes of 
solutions of the original system. This basic idea was developed and implemented recently 
by many authors, among others G. Bluman and J.D. Cole [10J, P. Olver and Ph. Rosenau 
[TTj . W. Fushchych [12], P. Clarkson and P. Winternitz [T3]. For a comprehensive review 
of this subject see the chapter by P. Olver and E.M. Vorobev in [T4], (Vol. 3, Chapter 
11). 

The method proposed in this work is distinguished by the specific choice of multiple 
differential constraints (DCs) of the first order, compatible with the initial system (11.61) 
for which the invariance criterion is identically satisfied. Consider a set of vectors 

£ a (u) = (u), . . . , e a (u)) T , a = 1, . . . ,p - 1 (2.1) 

which are orthogonal to a given wave vector A 

A M -e = 0, /i = l,...,p, o = l, ... (2.2) 

Note that £ a is not uniquely defined and the set {A, £i, ... , £ P -i} forms a basis in the 
space of independent variables E C R p . Multiplying equation f 1 1 . 3 2 j) by the vectors £ a and 
using the orthogonality property (12. 2p we obtain for any a = 1, ... , q and a = 1, ... ,p—l 



e(u(x))^ = 0. (2.3) 

This means that the functions u = (u 1 (x), . . . ,w 9 (x)) G U C R 9 are invariants of the 
vector fields (u(x))<9/<9a; M acting in E space. Hence the graph of the solution Y = 
|(x, u(x)j} is invariant under the vector fields 

X a = g(u)-jL, o = l,...,p-l (2.4) 

in the space of independent and dependent variables E x U C IR P x lR q . Note that the 
vector fields X a of the form (12.41) commute, i.e. they form an abelian distribution. 

Now, if the function u (x) is a solution of equations (ll.3ip . then u (x) is a solution of 
(11. 6j) if and only if the equation 

df a 

[A (/)I + Aj (/) A 1 (/)] -j- =0 (2.5) 

holds. This means that df /dr has to be an element of ker (A G (/)I + Aj (f)A % (/)). Equa- 
tions (12. 5p form an underdetermined system of first order ordinary differential equa- 
tions (ODEs) for /. Note that the differential constraints depend on the dimension of 
ker (Ao (/)I + A, {f)A l (/)). For example, if 

A (/)I + A, {f)A* (f) = , (2.6) 

then there is no differential constraint on the function / at all. 
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Thus, putting it all together, we can conclude that a function u (x) defined on a 
neighborhood of x = satisfies an equation of the form (ll.3ip 



u = / (a o (u)* + A-x) (2.7) 

for some / : R — > 1R 8 if and only if the graph T = { (x, u(x)) } is invariant under the 
vector fields X a of the form (12.41) . Such a function u (x) is a solution of the system (11.61) 
if and only if the ODEs (12.51) are satisfied. Note also that, due to the homogeneity of 
equations (12.51) . the rescaling of the wave vector A produces the same solution of the form 

flESU). 



Now we choose an appropriate system of coordinates on E x U space in order to rectify 
the vector fields (12.41) and find the invariance conditions which guarantee the existence of 
a rank-one solution d 1 . 3 1 [) of equation (jl.6p . We assume that one of the components of 
the wave vector A is different from zero, say Ai ^ 0. Then the independent vector fields 



X 2 = — - — — x = — - ^— (2 I 

dx 2 Ai dx 1 ' ' p dx p Ai dx l 



X 2 - > • • • > X P ~ 7^ > ( 2 - 10 ) 



have the form ( 12.4[) with the orthogonality property ( 12. 2D . If we change the independent 
and dependent variables as follows 

x 1 = r(x, u) , x 2 = x 2 , ... ,x p = x p , u 1 = u\ u 2 = u 2 , . . . , u q = u q , (2.9) 

then the vector fields (I2.8P take the rectified form 

d_ _d_ 

dx 21 '" ' p dxp 

and the corresponding invariant conditions are 

% 2 = 0, ... ,%p = 0. (2.11) 

So, in this new coordinate system on R p x R 9 , we subject the original system (11.61) to 
the invariance conditions (12.111) and produce an overdetermined quasilinear system of the 
form 

[A (u)I + A,(u)^(u)]u> =0, 
%2 = , . . . , u-xp = . 

with the general solution 

u(x) =f{x 1 ), (2.13) 
where the function / : R — * R 9 has to satisfy the ODEs (12.51) . 

Suppose now that the integral curve T of the vector field 7" (u)d/du a on U satisfies 
the ODEs (11.341) . Suppose also that the wave vector A (u) is pulled back to the curve 
T, i.e. /*(A). Then the functions A^ (u) become functions of the parameter r defined on 
the curve T. We denote /*(A At ) by A^ (r). The set (ll.24p and (ll.3ip of implicitly defined 
relations between the variables u a , x^ and r can be written as 

u a = f a (r), r = X fl (r)xr (2.14) 
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Expression (12.141) constitutes a rank-one solution (a Riemann wave) of system (11.61) and 
the scalar function r (x) is the Riemann invariant associated with the wave vector A (r). 

Now we present rank-one solutions of the MHD equations (1 1 . 1 [) - ( TT75T) which illustrate 
these theoretical considerations. Equations (11.121) . (I1.13P and (12.21) determine the char- 
acteristic directions A s , 7 S and £ a , each of which are of four types, and lead to different 
physical properties. We now discuss the rank-one solutions of MHD equations obtained 
from the reduced equations (12. 5p . The different symmetries associated with the vector 
fields X a of the form (12.41) lead us to four different types of solutions. The results can be 
summarized as follows. 



2.1 Simple entropic waves. 

— # 

We consider three types of entropic waves generated by the wave vector A. We denote by 
Ei, E 2 and E 3 the entropic waves related to wave vectors which are of one, two or three 
dimensions, respectively. 

1). For the entropic wave of type Ei the vector fields X a satisfying condition (12. 2p are 
given by 

{y.\)dt dx (v-X)dt dy (v ■ X) dt dz 

The solutions invariant under {X\, X 2 , X 3 } have the form 



\H 



2 



P = p(r), p(r) + —^-=p , v = v(r), H = H (r) , (2.16) 
where the Riemann invariant is 

r = A(r) -x- (A(r) • v)t , (2.17) 
p is an arbitrary constant, the vectors v and H have to satisfy 

'fxf\-ii = , (2.18) 
ar ar J 



and the wave vector A is given by 



t H dH . 
A = ^x— . 2.19 
\H\ 2 dr 



The relations (12.181) and fl 2 . 1 9 [) imply that the entropic wave solution (12.161) is constant 
on the planes perpendicular to A. Thus the entropic wave E\ is a plane wave which 
propagates in an incompressible fluid. The quantities p, p, v and H are conserved along 
the flow. The Lorentz force F m associated with the entropic wave E\ is given by 

F m = - l -V[\H\ 2 ], (2.20) 
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and is cancelled out by the hydrodynamic pressure, Vp, so that the fluid is force-free. 
Consequently, by virtue of Kelvin's theorem, the circulation 

v-d£ (2.21) 



c 

around a fluid element is preserved [13]. Note that F m -v ^ 0, which indicates the existence 
of coupling between hydrodynamic and magnetic effects. 

2). The entropic wave E 2 is characterized by (11.201) and the corresponding vector 
fields X a satisfying condition (12.21) are 

^ 1 d d ^ d X ^ (2 22) 

v dt v dx dy ' dz 

The solutions invariant under {Xi,X 2 } have the form 
1 
2 



p = p(r), p(r) + \\H\ 2 = p , v = u (r)e*i +v (r)e 2 + w(r)e 3 , H = H(r)e 3 , (2.23) 



with the constraint 

u (r) + v (r) = U , (2.24) 

where p Q and U Q are arbitrary constants, and p, p u, v, w and H are arbitrary functions 
of the Riemann invariant 

r = x + y -U t. (2.25) 
The Lorentz force related to the entropic wave E 2 is given by 

Fm = ~ l -V[\H\ 2 ] =Vp, F m -v = 0, (2.26) 

and therefore there is no coupling between hydrodynamic and magnetic effects. The 
entropic wave E 2 propagates in an incompressible fluid without vorticity. 

3). An entropic wave E3 is invariant under the vector field 

d d d d ,„ „. 

x = 7r+ +u 7T + v ir + w ir- 2 - 27 

ot ox ay oz 
Solving the MHD system (ll.ll) - (11.5j) we obtain 



\H 



2 



p = p(r), p(r) + ^-= Po , v = v(r), H = H (r) , (2.28) 

where the following constraints 

u (r) + v (r) + w(r) = C a , H x (r) + H 2 (r) + H 3 (r) = H Q , (2.29) 

^ = 0, V[p+ 1 -\H\ 2 ]=(H-V)H 

hold. Here C Q and 7i D are arbitrary constants. By taking into consideration the eigenvec- 
tor 7^3 = ^7 P , 0,0,C)V we obtain the solutions 

p = p(r), p = p , v = v , H = H Q , (2.30) 

where p Q is an arbitrary constant, v c and H are constant vectors and p is an arbitrary 
function of the Riemann invariant 

r = x + y + z -C t. (2.31) 

It is well known [15] that, after an appropriate Galilean transformation, the solution given 
by (I2.30p and (12.311) becomes stationary. 



9 



2.2 Simple Alfven waves. 



The vector fields X a associated with the Alfven wave A 6 are 

Xi = -^i— + — , X 2 = -^- + — , X 3 = -^- + — , (2.32) 
A D <9t <9x ' A D dt dy ' \ Q dt dz 1 

where 

A Q = e ^ H ' X) -v- A, e = ±l. (2.33) 
The solutions invariant under {Xi,X 2 ,X 3 } are given by 

P = Po, P = Po, v=^ + v , |#| 2 = W*, e = ±l, (2.34) 

where p Q , p Q , v D and 7i are arbitrary constants and H is an arbitrary vector function of 
the Riemann invariant 

r = A(r) -x+ (A(r) • v Q )t . (2.35) 
The vectors v and H satisfy the constraints 

dv -> dH -» 

— . A = 0, — -A = 0. 2.36 

ar dr 

The Lorentz force _F m associated with the Alfven A e wave is 

F m = (H-X)^. (2.37) 

The solution ( 12.341) . along with ( 12.371) . describes the well known properties of the Alfven 
waves [16], namely, they have no group velocity (the solution (12. 34ft is stationary, after 
an appropriate Galilean transformation) and propagate in an incompressible fluid. The 

Lorentz force F m acts on the fluid transversally to the direction of the wave vector, A, so 
that the magnetic field lines twist relatively to one another, but do not compress. Finally 
the fact that F m • v = indicates that there is no coupling between hydrodynamic and 
magnetic effects. 



2.3 Slow and fast simple magnetoacoustic waves. 



We consider now the slow S and fast F magnetoacoustic waves which are invariant under 
the vectors fields 

X l = -^— + — , X 2 = -^— + — , X 3 = -^— + — , (2.38) 
A D dt dx ' A D dt dy ' A D dt dz ' 

where 

\ = 5 s /f\\\ -v-\ s/f . (2.39) 

The functions Ss\X\ and Sp\\\ are given by equations (11.171) and f II . 1 8 [) . respectively. We 
denote them by Ss/f |A| and the wave vectors X s and X F by A 5//F . Under the conditions 
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(12.3§p the reduced system (12.51) takes the form of the following nonlinear system of ODEs 



dp 


rj(r)p 


dr 


S 2 


dv 


77(f) 


dr 


$S/F 


dH 




dr 


= i] (r) 



r) 2 

°S/F 



(H-X s / F f 
\X S / F \ 2 p 



\x s / F \ 

H-(H- X s / p ) 



x s / F 
\x s / F \< 



(2.40) 



where 

r/(r) = a" 



$S/F 



\H\ 
P 



, H- \X S ' F 



x 



dH 

dr 



= 0, H 



= 0, P=A oP K 1 



(2.41) 

and A Q is an arbitrary constant. The determination of the general solution of system 
(I2.40P is a difficult task. To simplify it we decouple the system (12.401) by considering two 
separate cases of the fixed orientation of the magnetic field H relatively to wave vectors 
X s and X F . 

1). The magnetic field H is orthogonal to X S ^ F . 

Here the slow magnetoacoustic wave S becomes an entropic wave E\ . The correspond- 
ing solution for the fast magnetoacoustic wave F is 



p = p ( r ) , p = A p K , H = pH a 
where p is an arbitrary function of the Riemann invariant 



r = A • x + JnAop- 1 + p\H \ 2 -v-X" )t. 



(2.42) 



(2.43) 



Here H Q is a constant vector perpendicular to the direction of wave propagation X F which 
is also constant. The velocity of the flow is parallel to X F and is expressed in terms of the 
density 

v = ev (p) X F , £ = ±1 , (2.44) 
where the norm of v is given by 



' \[A~ \\f(3 p + 1 — arc tanhi//3 p + 1]] for re = 1, 
v /2A^Tl) v /p for re = 2, 

(«-2)V^ 



V K A Q 



(re-1) 



Vp*' 1 + Pop 



Vp k ~ 2 + Po 



2 F 1 [a 1 ,bi,c 1 , 



do 



(2.45) 



for re ^ 1,2. 



The function 2-^1 (cli, bi, Ci, z) is a hypergeo metric function of the second kind of one 
variable z = —p K ~ 2 / /3 , where 

p ° = Tjt ' (2 ' 46) 
11 



with parameters a\ = 1/2(k — 2), b\ — |, Ci = 1 + 1/2(k — 2). Equation (I2.46P gives the 

ratio of the magnetic pressure ||-ff | 2 to the hydrodynamic pressure p ~ A a . When (3 Q is 

of the order of one or larger, the flow will be affected noticeably by the magnetic fluid H. 
If P <C 1, the opposite is true. 



The Lorentz force F m related to the fast magnetoacoustic wave F is 

1_ r , ,t,oi , ,o dp 



[2.47) 



and therefore the fast magnetoacoustic wave moving perpendicularly to H causes com- 
pressions and expansions of the distance between the lines of force without changing their 
direction. Since Vx v = 0, the fluid has no vorticity and the fact that F m is a conservative 
force implies, by virtue of Kelvin's Theorem, that the circulation (12.211) is constant. 



2). The magnetic field H is parallel to X S ^ F . 

For §a = | H | 2 / p < a, the fast magnetoacoustic wave F becomes simply an acoustic 
wave that propagates along a constant magnetic field H Q and there are no magnetic effects 
[TT] . The slow magnetoacoustic wave S (or fast F if 5 a > a) is more interesting: 

-| -1/(k+2) 

, p = A oP K , H = H X(r). (2.48) 



P 



Ik + 1 



(1+2/k) 



2/5 G J Po 

Here A Q and H Q are arbitrary constants and p is a function of the Riemann invariant 



r = A (r) ■ x + | A (r) | (5a — v sin 9)t . 
The flow velocity takes the form 

v = £v[sin 9 A (r) — cos 6 Aj_(r)] , e = ±1 



(2.49) 



(2.50) 



where the index _L specifies the direction perpendicular to A and the angle 9 = Z(A, 7) is 
defined by 



tan 9 



The norm of v is given in terms of p by 



SI 



-Pop- K - 



(p) = 2H oP 



P lK + Po - r, nrs 2F1 ( a u 61, ci ; 



1 + 2k 



2k 



(2.51) 



(2.52) 



where the function 2 -Pi (ai,6i,ci;2) is a hypergeometric function of the second kind of 
one variable z = —p 2K /f3 2 , where (3 a = H 2 /kA q , with parameters a\ = —(1 + 2k)/4k, 
61 = |,ci = 1 - (l + 2«)/4«. 

For /3 ~ ^i/ a2 ~ 1) when the magnetic pressure and hydrodynamic pressure are 
comparable, the magnetic field lines undergo compressions and expansions resulting from 
the behaviour of the longitudinal component of 7 along the wave vector A (similar to 
acoustic waves in hydrodynamics). In addition the tensions produced by the Lorentz 
force F m are of the form 

F m = ^-X±(r). (2.53) 
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The slow magnetoacoustic wave S is often called the "compressional Alfven wave" [IB] . 
If p o 3> 1, the magnetic field is so strong that 

(A-7) = ^^0 ^ 0--, (2.54) 

and we obtain the simple Alfven wave A £ described in Section 2.2 which propagates in 
an incompressible fluid. This wave is called the "torsional Alfven wave" [IB] . 

The simple wave solutions of the MHD equations described above are not new. They 
have been obtained in the past, mostly by the generalized method of characteristics (see 
e.g. [6]). It is worth mentioning that the treatment of MHD equations, or any hyperbolic 
equations for that matter, by the classical symmetry reduction method has never yielded 
simple wave type solutions. It has only become possible through the use of the conditional 
symmetry method, the variant of which we present here. Interestingly it delivers the 
solutions in a closed form, lending themselves easily to physical interpretation. The GMC, 
on other hand, usually provides solutions in much more complex form and a good deal 
of guesswork is often involved in bringing them to more serviceable form. Thus, it seems 
that there exists a certain affinity between conditional symmetries and the structure of 
simple wave solutions of hyperbolic equations. This fact is of more practical consequence 
when we proceed to multiple wave solutions. 



3 Conditional symmetries and double waves in MHD equations. 

Now we generalize the above construction to the case of superposition of many simple 
waves described in terms of Riemann invariants. For this purpose we fix k linearly inde- 
pendent wave vectors A 1 , ... , X k , 1 < k < p, which correspond to the following Riemann 
invariants 

r a (x,u) = A*(u)x", s = l,...,k. (3.1) 

The equation 

u = /(r 1 (x,u), ... ,r fc (x,u)) (3.2) 

defines a unique function u (x) on a neighborhood of x = and the Jacobian matrix is 
given by 




with l,j — 1, . . . , k ; yu — 1, . . . ,p ; a — 1, . . . , q where the matrix 

is assumed to be invertible. This assumption excludes the gradient catastrophe phe- 
nomenon for the function u. Note that the rank of the Jacobian matrix (13. 3p is at most 
equal to k. If the set of vectors 

£ (u) = (£(u), ... ,^(u)) T , = 1, ... ,p-k (3.5) 
satisfies the orthogonality condition 

A*-e = 0, s = l,...,k, a = l,...,p-k (3.6) 
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then we get 

eKx))^ = 0. (3.7) 
This means that u 1 (x), ... , u q (x) are invariants of the vector fields 

e(u(x))^, a = l, ...,p-k (3.8) 

in the space of independent variables E C R p . Hence the graph of the solution (13.21) 
T = { (x, u(x)) } is invariant under the vector fields 

X a = e(u)^, a = l,...,p-k (3.9) 

acting on the space of independent and dependent variables E x U C R p x R 9 . 

If u (x) is a g-component function defined on a neighborhood of x = such that the 
graph T = { (x, u(x))} is invariant under all vector fields (I3.9P with the property (13.61) . 
then u (x) is the solution of (13.21) for some function /. The functions r 1 , . . . , r h , u 1 . . . , u q 
constitute a complete set of invariants of the abelian algebra generated by the vector fields 
(13. 9p . Substituting expression (13.21) into system (II. 6p we obtain the reduced system 

r 1 W)^^(u(x))^(u(x))|^(r(x,u(x))) =0. (3.10) 

Under these circumstances we can state the following. 

Proposition. A function u(x) defined on a neighborhood of x = satisfies an equation 
(13. 2p for some function / : R fc — > R 9 if and only if the graph V = j(x, u(x)) } is invariant 
under all vector fields (13.91) . Such a function u(x) is a solution of the system ( 11 .61) if and 
only if the partial differential equations (13.101) are satisfied. 

Now we show that a proper change of variables allows us to rectify the vector fields 
(13. 9p and, consequently, to derive a reduced system which admits multiple wave solutions. 
If the k by k matrix 

A=(A]), l<i,j<k (3.11) 
is invertible, then the independent vector fields 

Xk+1 = " ^ (A ~^ Xk+1 ^ ' • • • ' Xp = ^ ~ £ (A ~^ A ^ (3 - 12) 

i,j=i i,j=i 

have the required form (13. 9p for which the orthogonality conditions (13.21) are satisfied. 
The change of independent and dependent variables 

x 1 = r x (x,u), ... ,x k = r k (x, u), x k+l = x k+1 , . . . ,x p = x p , u 1 = u 1 , . . . , vfl = u q 
allows us to rectify the vector fields X a and we get 



d d 

x ™=dP+i>-> x ' = mi- (3 - 14) 
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The corresponding invariance conditions are 

%fc+i = , . . . , u-xp = . (3.15) 
The general solution of the invariance conditions H3. 15f) is given by 

u(x) = /(x\ ...,x fe ), (3.16) 

where / : IR fc — > R 9 is arbitrary. The system ( 11. 6ft is subjected to the invariance conditions 
( 13.15!) and, when written in terms of new coordinates (x, u), takes the form 




(3.17) 



Note that for k > 2 the system (I3.17P is more sophisticated than its analogue f|2.12j) . This 
is due to the presence of the k by k matrix $. 

As it was shown in [1] (page 885, Theorem 2), that the symmetry criterion for the 
existence of group invariant solutions of the form (I3.16P (z.e.rank-& solutions) for the 
overdetermined system (13.171) states that 

pr (1) X a A = 0, a = l,...,p-k (3.18) 

whenever equations A = are satisfied. 

We complete our construction of rank-A; solutions of system (II. 6p by pulling back the 
wave vectors A 1 , ... , A fe to the k- dimensional submanifold S C U obtained by solving 
the system ( I3.17p . Then the wave vectors A 1 (u), ... , A fc (u) become functions of the 
parameters r 1 , . . . , r k . The set ( 13. ip and (13.21) of implicitly defined relations between the 
variables u a , x M and r 1 , ... ,r k can be written as follows 

u = /(r), r s = A*(r)^, s — 1, ... ,k (3.19) 

where r = fr 1 , . . . ,r fe ). Note that the rank of function u (x) is at most equal to k. In 
particular solutions (13 19j) include Riemann k- waves (as presented in Section 4). 

The construction procedure outlined above has been applied by us to the MHD equa- 
tions (ll.ll) - (ll.5p in order to obtain rank-2 solutions representing nonlinear superpositions 
of two simple waves. The results of our analysis are summarized in Table 1, which shows 
the possibility of existence of these solutions obtained from different combinations of the 
vector fields X a . 

The remaining part of this section is devoted to the discussion of some more interesting 
solutions. 

We denote by EiEj, A e A e , A e Ei, FF, FEi, . . . , i, j = 1, 2, 3, the solutions which are 
the result of nonlinear superpositions of given waves. The indices i,j are related to the 
dimension generated by the wave vectors A s . Moreover we denote by r and s the Riemann 
invariants (which coincide with group invariants of X a ) of the waves under consideration. 
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3.1 Double entropic waves. 



Analyzing the double entropic waves E^E^ we consider separately two cases. 

1) . In the first case two wave vectors A 1 and A 2 are located on the plane spanned by 
the unit vectors e*i and e.2- Hence we have 

A 1 = ( - (A 1 • v), A 1 ), A 2 = ( - (A 2 • v), A 2 ) , (3.20) 

A 1 = (cos<^, siny), 0) , A 2 = (cos^, sin6>, 0) . (3-21) 

Here <p and 9 are functions of the Riemann invariants s and r which are to be determined. 
The general solution of the type E\E\ invariant under the vector fields (13. 9p 

v at v ox Oy oz 

is determined by the following system of equations 
j i 

— {p,p,v,H} = 0, V-v = 0, Vp = --V|#| 2 + (H ■ V)#, (#-V)v = 0, V-# = 0, 

(JjL £ 

(3.23) 

where d/dt = d/dt + (v • V); p, p, v and H are arbitrary functions of the Riemann 
invariants 

s = X 1 (s,r) ■ (x- vt) , r = A 2 (s,r) • (x- vt) . (3.24) 

The reduced system f !3.23|) describes the propagation of a double entropic wave E^E^ in an 
incompressible fluid for which the quantities p, p, v and H are conserved along the flow. 
The fluid is force free since we have an equilibrium between the gradient of hydrodynamic 
pressure and the Lorentz force F m which is the sum of the gradient of magnetic pressure 
\H\ 2 /2 and a tension force (H-V)!!. The constraint (H-'V)v = implies that the velocity 
of the fluid v remains constant along H. The double entropic wave produces vorticity and 
the circulation f)2.2ip of the fluid is preserved since dv/dt = 0. If we impose the condition 
that the Lorentz force depends only on the magnetic pressure, then we obtain the solution 

p = p(s,r), p(s,r) + ^\H(s,r)\ 2 = p , v = w{s,r)e 3 , H = H(s,r)e 3 , (3.25) 

where p, w and H are arbitrary functions of the Riemann invariants s = A 1 (s, r) ■ x and 
r = A 2 (s,r) • x; p is an arbitrary constant. The double entropic wave propagates in an 
incompressible fluid in which the flow is one-dimensional and stationary without vorticity. 

2) . In the second case we fix the direction of the wave vector A 1 along e 3 and we let 
the wave vector A 2 circulate on the plane perpendicular to A 1 . The entropic wave vectors 
have the form 

A 1 = (-«,, 0,0,1), A 2 = ( - (A 2 ■ v), A 2 ) , (3.26) 

where A 2 = (cos#, sin 9, 0) and 9 is a function of the Riemann invariants s and r which 
are to be determined. The corresponding vector fields, (13. 9p . are given by 

X 1 = -tan^— + — , X 2 = — + (u + vt&n9)— + w— . (3.27) 
ox oy at ox Oz 
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We can now consider two types of MHD solutions, invariant under {Xi, X2}, which depend 
on the direction of the magnetic field H with respect to the wave vectors A 1 and A 2 . 



—> — » . — * 

2. a). The magnetic field H is perpendicular to the wave vectors A and A . The 
solution is given by 

P = p{s,r), p{s,r) + ^\H(s,r)\ 2 = p a , v = v(s, r) Aj_(r) + w(s) e 3 , 
H = H(s,r)X 2 ± (r), (3.28) 

where p a is an arbitrary constant; p, v, w and H are arbitrary functions of their arguments, 

the unit vector A^ = (— sin 9 (r), cos 9 (r), 0) is perpendicular to the wave vector A 2 and 
9 is an arbitrary function of r. The Riemann invariants s and r are given by 

s = A 2 (r)-x, r = z-w(s)t. (3.29) 



2.b). The magnetic field H is perpendicular to the constant wave vector A 2 . Here the 
solution takes the form 

p = p(s,r), p(s) + ^\H{s)\ 2 = p , v = v(s)e + w{s)e 3 , H = H±(s)e + H 3 (s) e 3 , 

(3.30) 

where p Q is an arbitrary constant, the constant vector e Q = (— sin 9 a , cos 9 a , 0) is per- 
pendicular to the wave vector A 2 , p is an arbitrary function of two Riemann invariants 
s and r; and v, w, H± and H 3 are arbitrary functions of the Riemann invariant s. The 
invariants s and r have the form 

s = A 2 -x, r = z — w(s)t . (3.31) 

For both cases 2. a) and 2.b) the Lorentz force is cancelled by the gradient of the hy- 
drodynamic pressure. Thus the fluid is force free. By virtue of Kelvin's Theorem the 
circulation (12.211) of the fluid is preserved. 



3.2 Double Alfven wave AA. 



We discuss the superposition of two Alfven waves A 6 for which the wave vectors have the 
form flL22D 

A 1 = ( Aq , A}, A 2 , A3 ) , A 2 = ( A 2 , A 2 , A 2 , A 2 ) , (3.32) 

where 

xi (H- ■A') i = l,2, e = ±l. (3.33) 

We assume that the wave vectors A 1 and A 2 are linearly independent. The corresponding 
vector fields (13. 9p are given by 

[\\\\ - A*A 2 ) d - A*Af ) d d_ 

- dt + (XlXl - A^A 2 ) dx + Oy ' [6 - 6 } 

(A}Aj - A 3 A 2 ) d (A3A 2 - A^A 2 ) d d_ 



(A^A 2 - X\Xl) dt (XlXj - X\Xl) dx dz 
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The solutions invariant under {Xi,X 2 } are 

P = Po, P = Po, v(s,r) = -^=H(s,r) , \H\ 2 = V 2 , e = ±1 , (3.35) 



where p D , p Q and are arbitrary constants, v and if are vector functions of the Riemann 
invariants 

s = X 1 (s, r)-x, r = A 2 (s, r)-x. (3.36) 
The magnetic field H must satisfy the constraints 

— -A^O, — -A 2 = 0. (3.37) 
as dr 

The double Alfven wave AA propagates in an incompressible and stationary fluid which 
flows along the lines of force of the magnetic field H . 



3.3 Double Alfven entropic waves AEx- 

We now consider the superposition of the Alfven wave A 6 with the entropic wave Ex- The 
wave vectors (I1.19P and (I1.22I) are given by 

A 1 = ( A* , A}, A 2 , A3 ) , A 2 = ( A 2 , A 2 , A 2 , A 2 ) , (3.38) 

where the vectors A 1 and A 2 correspond to the Alfven X A and the entropic X E wave vectors, 
respectively, and 

A^ =e^-^l-v-A 1 , e = ±l, A 2 = -A 2 -v. (3.39) 
\[P 

The invariant solutions under the vector fields (13. 9p . 

(AjA| - AgA 2 ) d (XlXl - XjX 2 ) d d_ 

1 " (AJA? - X\Xl) dt + (ASA? - AJA?) &c + dy ' l ' j 

(AjAf - A*A 2 ) g (A3A 2 - A^Aj) d d 

2 ^ (KXl - X\Xl) dt + (KXl - AJA*) + dz ' 

take the form 

P = p(r) + ^ 2 (r)= Po , |tf| 2 = H 2 (r) , 

-» -» 
+ <p(r), # = a(s,r)^+/3(s,r)(^ + ^(r) , e = ±1 . (3.41) 



They are expressed in terms of the Riemann invariants 

s = X 1 {s,r)-x, r = X 2 (s,r) ■ x- (A 2 (s, r) ■ v)t , (3.42) 

where p Q is an arbitrary constant, p, p and TC are arbitrary functions of their arguments. 
Furthermore the vector functions <p and ip must satisfy the algebraic relations 

0-X 1 = O, <pT-A fc = 0, fi-X k = 0, 0-X 2 = O, $-X 2 = 0, ^-A 2 = 0, (3.43) 



18 



for k = 1, 2. In expressions (13.411) and (I3.43[) we have used the following notation 



<P 



<ir 2 



<ir 3 



dr 



(3.44) 



The function a is expressed in terms of the function 3, the vector functions (p and if) and 
their derivatives 

-[(^-^)/3 + (^-VT)]± VA 



a 



(3.45) 



where the solution (13.451) is real if the discriminant A is positive, i.e. 



A 



(<P-<P)P +(<?•'$) 



m 



U\r) 



> 0. 



(3.46) 



To calculate the function (3 we must solve the following nonlinear ODE 

.j]Q3 r • •■ --i r • -> - • ^ •■ i 

x ip) • ^-^7+ (0 x £)• £ 2 + x # £ +(£ x <P) ■ i> ~\0\' 2 {0- <P)[(<P x ■ 0] P 



0\- 2 [{0x^)-ip] (^-^tVA +[(^x^.fl = ) (3.47) 



where the coefficients are expressed in terms of <p, if) and their derivatives. Analyzing the 
solutions (I3.4ip . we notice the synergetic effects which result from the superposition of 
the two simple waves. This superposition constitutes the double wave AE\. We have an 
equilibrium between hydrodynamic and magnetic pressures resulting from the entropic 
wave E\. In addition we observe that the density p and the pressure p are unchanged 
by the Alfven wave A 6 . Both waves modify the flow direction that is not parallel to the 
magnetic field H, in which case the flow is incompressible and rotational (i.e.Vxv^O). 
However, the circulation of the fluid is not conserved because of the existence of tension 
forces originating from the Lorentz force (due to the contribution of the Alfven wave A e ). 



3.4 Double magnetoacoustic waves FF. 



In the study of the double magnetoacoustic waves FF we limit our analysis to the situation 
in which the magnetoacoustic wave vectors A 1 and A 2 are orthogonal to the magnetic field 
H. We consider two cases. 



1). We assume that H = (0, 0, if) and the wave vectors are of the form 



A 1 = (6 F -u, 1,0,0) 



A 



(5 F -v,0, 1,0), 



where 



up \H\ 

h J — - , e = ±1 . 

P P 



Then the corresponding vector fields X a are given by 



X 2 



d_ 

dz 



(3.48) 



(3.49) 
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The solutions invariant under {Xi,X 2 } are 



1 



P = p(s,r), p = A Q p K , v= v(s,r) + -[f(r)+g(s)] 



H 



pH e 3 



e*i + v(s, r)e 2 + w(s - r)e 3 , 

(3.50) 



where A D and H Q are arbitrary constants; p, v and w are arbitrary functions of their 
arguments. The Riemann invariants are of the form 

s = x + (s F -^[f(r)+g(s)] -vjt, r = y + {5 F -v)t, (3.51) 

and the functions / and g are given by 

for k = 1, 



f(r)-g[s)=e< 



WA 



\/ (3 Q p + 1 - arc tanh [\J (3 Q p + 1] 
±yJlAo{p o + l)p for k = 2, 



1 + 



2-Pi ai,6i,ci, 



-P 



K-2' 



(3.52) 



for k ^ 1,2 



where e = ±1, the function 2 F\ (a\, & a , c\\ z) is a hypergeometric function of the second 
kind of one variable z = —p K ~ 2 //3 and (3 Q = H^/k,A , with parameters a± = 1/2(k — 2), 
6i = i,ci = 1 + 1/2(^-2). 



The Lorentz force associated with the double wave FF is 
F m = -\V[\H?]=-Hlp 

Consequently the double wave FF causes compressions and expansions of magnetic field 
lines in two directions e% and e 2 , perpendicular to the magnetic field H. Since F m is 
a conservative force (i.e. it can be derived from the gradient of the magnetic pressure 
||i/| 2 ), the circulation of the flow (12.211) is preserved. 



dp _ dp _ 
as or 



(3.53) 



2). We consider the one-dimensional case of a superposition of two magnetoacoustic 
fast waves which propagate with local velocities = ^ + eu, e = ±1. So we have 



A 1 = (8 F -u, 1,0,0), A 2 = (-(5 F +u), 1,0,0) 



where 



up 
P 



The vector fields, f)3.9p . take the forms 



1 dy ' 



\H\ 



oz 



(3.54) 
(3.55) 

(3.56) 



Here we assume that the unknown functions u = (p,p,v,H) are some functions of the 
Riemann invariants s = x + (5 F — u)t and r = x — (5 F + u)i. Substituting these functions 
into the MHD equations fll.ip - fll.5l) we obtain the solutions 



P 
H 



A Q p K , v = ^[/(r) +g(s)] e x + v(s + r)e 2 + w(s + r) e 3 



p\s,r), p 
pH [ cos ip(s + r)e 2 + sin <^(s + r) e 3 ] 



(3.57) 
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where A D and H Q are arbitrary constants; p, v , 10 and are arbitrary functions of their 
arguments. The Riemann invariants are of the form 



x 



+ (*F-i[/(r)+0(*)])t, r = o?- (fc + i[/(r)+0(a)])t, (3.58) 
where the functions / and g retain the form (13.521) . 

For the particular 2 the solutions (13.571) can be expressed as 

[f(r)-g(s)] 2 . l [f/u Ml . ft „ - 

P = 16 ^ 2A +H2 ^ , P = A p, v = - [f(r) + g(s)\ e x , H = pH e 3 , 

where A and H Q are arbitrary constants, / and g are arbitrary functions of the Riemann 
invariants r and s, respectively. These invariants are given by 

s = x-^{3g + f)t, r = x-~(3f + g)t. (3.59) 



The Lorentz force is 

I _Tr)n fin 

e x . (3.60) 



F m = - l -V[\H\ 2 } =-H 2 oP 



dp dp 
ds dr 



Thus the double wave FF propagates in the fluid as the combination of a compressional 
wave associated with the arbitrary function fix — (§f + u)t) and an expansion wave 
associated with the arbitrary function g(x + (5f — u)t). Since V x v = 0, the fluid is 
irrotational. This particular case (k = 2) provides an example of the MHD generator 
principle, i.e. conversion of mechanical energy into electromagnetic energy. The fluid 
flows with a velocity v = ue\ and the magnetic field H lies along the direction eV This 
results in an electric current J which takes the form 



dg_ _df_ 

ds dr 



e 2 . (3.61) 



Finally, if the double wave FF propagates along the magnetic field H, we obtain a 
double acoustic wave in hydrodynamics |17j . 



3.5 Double magnetoacoustic entropic waves FE%. 

We assume that the wave vector \ F of the fast magnetoacoustic wave F and the wave 
vector X El of the entropic wave E% are orthogonal to the magnetic field H. We study the 
following two cases. 

1). Consider the one- dimensional case of a superposition of an entropic wave E\ and 
a magnetoacoustic wave F propagating in opposite directions that are perpendicular to 
H = (Hi, H 2 , 0). We assume that their phase velocities are Af = — w and Af = <5p — w, 
respectively. We have 

\ F = (S F -w, 0,0,1) , X El = (-w, 0,0,1), (3.62) 
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where 



Sp = £ 

The vector fields, (j3.9p . take the forms 



np \H\ 2 
P P 



d_ 

dx ' 



£ = ±1. 



5 



9y 



(3.63) 



(3.64) 



We assume that the unknown functions u = (p, p, v, .ff) are some functions of the Riemann 
invariants s = z + (5p — u>)£ and r = z — wt. Substituting these functions into MHD 
equations fll.ip - fll.5D we obtain 



P = p(s) , p = A Q p K , v = a(r) x e 3 + to(p) e 3 
if = H Q p(s) [cos <p(r) e*i + sin <p(r) e 2 ] , 



(3.65) 



where t4 d and if are arbitrary constants; p, <p and a are arbitrary functions of their 
arguments. The Riemann invariants s and r are given by 



s = z- (w(p) - e [kAoP^-^ + F G 2 p] 2 ) t , r = z- w(p)t , e = ±1 

The function u> has the form 

/ A^[v / /^oP + 1 — arctanhfv^p + T]1 for k—1, 



(3.66) 



w 



{p) = 2{ yf^Ao 



(K-l) 



VP"" 1 + PoP 



(k-2)V^ 



2 Fi(ai,6i,ci, 



for^l,2. (3 - 67) 



The Lorentz force related to the double wave FEi is 

F m = -\nH?]=-Hl P f s e, 



(3.68) 



and represents the contribution of the fast magnetoacoustic wave F. The gradient of 
hydrodynamic pressure takes the form 



Vp=kA oP ^^- 



ds 



e 3 



(3.69) 



It follows that the magnetic field lines undergo compressions and expansions along the 
direction of the wave vector X F . On the other hand the entropic wave E\ contributes to 
the flow vorticity 

Vxv = e 3 xfce 3 j. (3.70) 
Since F m is a conservative force, the circulation (I2.2ip of fluid is preserved. 



In the special case k — 2 the solution (13.651) is reduced to the following expressions 

v = d(r) x e 3 + 2ey/C 2 p(s) e 3 , e = ±1 , 



P = p(s) , P = A(r)p 2 , 
H = p(s) TC(r) [cos <p(r) e± + sin <p(r) e 2 \ 



(3.71) 



where C 2 in an arbitrary constant. The arbitrary functions A(r) and 1~i{r) satisfy the 
algebraic relation 

2A{r)+H 2 {r) = C 2 . (3.72) 
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Thus the entropic wave E\ has an influence on the Lorentz force, but does not affect the 
equilibrium between the hydrodynamic and magnetic pressures. 

2). We now consider the case for which both wave vectors \ F and X El are orthogonal 
to the magnetic field H = (0, 0, if) 



{8 F -u, 1,0,0), A £l = (-v, 0,1,0) 



where 



np H 2 
.9 P . 



In this case the vector fields, ( 13. 9p . take the forms 

_ 1 9 (u-8 F ) d d 
1 v dt v dx dy 



±1 . 



dz ' 



(3.73) 
(3.74) 

(3.75) 



The solutions associated with the vector fields, (13.751) . exist and for k = 2 they are given 
by 



p = p(s), p = A(r)p 2 , v= 6(r) - 2ey/C 2 p(s) e\ + v a e 2 + w(r)e 3 , H = p(s)H(r)e 3 , 

(3.76) 

where e = ±1, C 2 and t> are arbitrary constants; p, b and w are arbitrary functions of 
their arguments. The Riemann invariants are 



x 



+ \ 3£^ C 2 p(s) - 6(r)J t , r = y-v Q t, 

Hence the functions A(r) and 7Y(r) satisfy the algebraic relation 

2A(r)+H 2 (r) = C 2 . 

The Lorentz force takes the form 



±1 



(3.77) 



(3.78) 



(3.79) 



which is the result of both waves' contribution, but the entropic wave E% does not disturb 
the equilibrium between hydrodynamic and magnetic pressures (see ( I3.76P and (I3.78P ). 
The compressions and expansions of the fluid must take place in the direction of the fast 
magnetoacoustic wave propagation X F , otherwise the flow would be incompressible along 
the directions e 2 and e%. The fluid vorticity is given by 



dw _ db _ 
Vxv = — ex - — e 3 . 
dr dr 



(3.80) 



This indicates that the vorticity depends only on the entropic waves E\. Kelvin's Theorem 
ensures that the circulation (12.211) of the fluid is preserved. 



4 Concluding remarks. 



In this paper we develop the conditional symmetry method as a tool for recovering the 
solutions of hyperbolic systems in the form of Riemann waves and their superpositions. 



23 



This idea was firstly outlined in pQ where, among others, the symmetry criterion for 
the existence of group invariant rank-fc solutions was derived. The progress made here 
consists in adapting the conditional symmetry method procedure to the solutions based 
on Riemann invariants. Its most important element is the introduction of the requirement 
that the graph of the solution be invariant under vector fields, (13.90 , with the orthogonality 
property (13.61) . When this takes place, the appropriate change of variable, (13.131) . allows 
us to rectify the vector fields (13. 9p and, consequently, determine the invariance conditions 
(I3.15P which are to be imposed on the initial system of equations (II. 6p . 

Riemann waves and their superpositions have been studied so far only in the context 
of the generalized method of characteristics [19] . This approach requires imposing certain 
conditions on vector fields 7 S and X s of the waves entering into an interaction. In order 
for us to comment on its relation to the conditional symmetry method we outline here 
the basic assumptions of GMC 

A form of solution, called a Riemann fc-wave, is postulated for which the matrix of the 
tangent mapping du is given by 

^(x) = £nx)7 s Q (u)A*(u), (4.1) 

s=l 

where £ s 7^ are treated as arbitrary functions of x and we assume that the vector fields 
71, ... , 7^ are linearly independent. We assume also that commutators of all vector fields 
7z and 7 S are linear combinations of these fields 

[lh Is] e span {7,, 7 S } , V I ^ s = 1, . . . , k . (4.2) 

If these conditions are satisfied, then, due to the homogeneity of wave relation (11.121) . we 
may change the lengths of the vectors ji and j s , so that the commutators of these vector 
fields vanish, i.e. 

[71,7.] = 0, VZ^s = l, ... ,k. (4.3) 

The vector fields 71, ... ,7^ form an abelian distribution on the space U. There exists a 
parametrization of the integral surface S in U tangent to these fields, 

S : u = /(r\ ... ,r fc ), (4.4) 

such that 

J^ = 7s, Vs^ = l, ... ,k. (4.5) 

The wave vectors X s (u) become functions of the parameters r 1 , . . . , r k . Consequently the 
differential of (14. 4p gives 

s=l A*=l 

which, together with the assumption (14.11) . leads to a system of exterior forms 

dr s = £ s (x)A* (r\ ... ,r k ) dx» . (4.7) 

It has been shown ([20], [21]) that the system (14.71) has solutions if the following conditions 
are satisfied 

dX s 

-— y G span {X s , X 1 } , Va^Z = l, ...,k. (4.8) 
or 1 
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Conditions ( 14. 3D and ( 14. 8 j) ensure that the set of solutions of the system (II. 6ft . subjected 
to (14. ip . depends on k arbitrary functions of one variable. It has been proved [21] that 
all solutions, i.e. the general integral, of the system (14 ,7p . under conditions ( 14. 8p . can be 
obtained by solving, with respect to the variables r 1 , ... ,r k , the system in implicit form 

A'(r 1 ) ...,r*)a* = ^(r 1 ,...,r*) ) (4.9) 

where ip s are arbitrary functionally independent differentiable functions of k variables 
r 1 , . . . ,r k . The solutions of (14 .9p are constant on (p — /c)-dimensional hyperplanes per- 
pendicular to wave vectors X s . 

As we can see, both methods discussed here exploit the invariance properties of the 
original system of equations. In the GMC they have the purely geometric character (14.41) . 
whereas in the case of CSM we make use of the symmetry group properties (13. 9p . The 
two approaches describe two different facets of the same geometric object. 

There are, however, two basic differences between the generalized method of charac- 
teristics and our approach. Riemann multiple waves defined by (14.11) constitute a more 
limited class of solutions than the rank-A; solutions postulated by the conditional sym- 
metry method. This difference results from the fact that the scalar functions £ s in (14.11) 
(describing the profiles of simple waves) are substituted in our case (in expression (13.31) ) 
with a k by k matrix $ which allows for much broader range of initial data. 

The second difference consists in the fact that the restrictions (14 .3p and (14.81) on the 
vector fields 7 S and X s , ensuring the solvability of the problem by the generalized method 
of characteristics, are not necessary in the conditional symmetry method. This makes 
possible for us to consider more complex configurations of simple waves entering into an 
interaction. 

In order to verify the efficiency of our approach we have used it for constructing rank-2 
solutions of MHD equations. The results obtained in this work coincide to a great degree 
with the ones obtained earlier by the means of the generalized method of characteristics. 
All our rank-2 solutions are in fact Riemann double waves (in the sense of definition 
(14. ip ). This is probably largely due to the fact that, to facilitate computations, we have 
introduced additional simplifications, assuming some specific configurations of the wave 
vectors X s included in superpositions. Most of the solutions obtained here are already 
known [§]. Some, however, namely solutions (15301) - (13321 . (15371) - (15351) and (13351) - (I53TD 
are new. It is worth mentioning that, at least in the case of MHD equations, our method 
proved to be much easier to implement than GMC. In particular the invariance condition 
(I3.15P provides us with a set of coordinates very convenient for integration of the reduced 
system ( I3.17p . 

Finally it has to be said that for many physical systems, e.g. nonlinear field equations, 
fluid membrane equations, there have been very few, if any, known examples of multiple- 
wave solutions. The version of the conditional symmetry method proposed here offers a 
new, and it seems promising, way to deal with this problem. 
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Table 1: Double Riemann wave solutions described by MHD equations, + indicates that double 
wave solutions exist; — indicates that double wave solutions do not exist, e = ±1. 



Waves 
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A^ 
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+ 


+ 


+ 




A* 
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+ 






F L 
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+ 




S b 
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